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ABSTRACT 


This thesis presents a method of reducing the time re- 
quired to accomplish ship design feasibility studies by 
coupling a naval surface combatant synthesis model with a 
general purpose nonlinear optimizer. Brief descriptions of 
optimization techniques and synthesis models are presented. 
The selection of design variables, constraints, and objec- 
tive function 1s presented via a design example. Further 
examples are given illustrating the ability to start with 
infeasible designs and proceed, through optimization, to 
feasible designs. Examples of ships developed using 
different objective functions are also presented. The 
synthesis/optimizer system provides results in a rapid and 
descriptive manner which compare favorably with existing 
designs and provides the naval architect with an efficient 
tool to use in the conceptual design phase of naval surface 


combatants. 
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I. INTRODUCTION 


A. BACKGROUND 

In today's environment of rapidly increasing costs, 
technological complexity, and growing threats, we must 
actively pursue ways in which to improve the effectiveness 
with which we apply the limited resources available to the 
design of naval combatants. Currently, the design process, 
from feasibility studies through keel laying of the lead 
Ship takes almost a decade to complete. There is growing 
pressure to accelerate this process wherever vossible in 
the design sequence. Inflation appears to be the primary 
motivator for this pressure; as we can look to painfully 
expenSive examples of construction projects, which when 
delayed for whatever reason resulted in increased costs 
without significant platform improvements as compensation. 

The naval architect and naval engineer involved in the 
design process today can turn to the high speed digital 
computer for a viable means of reducing the time required 
for the design process and thereby achieve significant 
monetary gains while maintaining or improving the quality 
of the design product. The goal of "best ship at least 
cost'' and the "computer revolution" have resulted in a 
very successful software development effort within naval 


research activities. Reference [1] provides a catalog of 
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computer-aided design and construction (CASDAC) programs 
currently used by the Naval Sea Systems Command. The pro- 
grams range from routines used for data management to 
complex synthesis model and arrangements packages. 

The design of naval vessels can be categorized into 
three major phases. These are: 

1. Conceptual Design 

2. Preliminary Design 

ee Contract Design 
At each phase, the ship is defined in greater detail than 
the previous phase. 

Conceptual design can be further subdivided into two 
components; 1) the feasibility study and, 2) the concept 
design. These too, reer levels of design detail with 
feasibility studies being the first and the crudest esti- 
mate of the ship, and concept design being the development 
and optimization of a single or several ships selected 
from feasibility studies. 

Mills, in [Ref. 2], defines feasibility studies as 
",..an estimate of the ship system level physical 
characteristics and cost related data for a design 
which represents a feasible solution to a specific 
set of performance requirements." 

Restated, a feasibility study is a shortcut estimate of a 
Ship's principal characteristics, machinery systems and 
various coefficients which when taken as a whole, are a 


feasible solution to the owner's specifications or desires. 
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Feasibility studies can be used for trade-off studies where 
new designs are evaluated; for subsystems trade-off studies, 
i.e. how different propulsions systems might work in a 
DD-963 hull, for evaluating changes in design standards and 
practices, or for determining an optimum feasible design 
with which to start the preliminary design phase. 

It is evident from the brief descriptions above that 
large numbers of feasibility studies are required to satisfy 
all of the designers needs. This demand for volumes of 
output and information has resulted in feasibility studies 
being successfully computerized in the form of synthesis 
models. A synthesis model is an engineering design 
procedure for converting a set of requirements into the 
moysical description of a ship which can satisfy those 
requirements. The ultimate result of using the synthesis 
model as a design tool is the ability to produce a far more 
detailed and accurate design earlier in the design sequence, 
thereby saving precious time and money and providing more 


reliable guidance in the design selection process. 


B. PROBLEM STATEMENT 

The main problem, simply stated, is that the design of 
naval surface combatants must be accomplished as quickly and 
as inexpensively as possible to prevent inflationary over- 
runs and technology lag. The past three decades have 


provided numerous examples of the consequences of "too long 
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miadetoo much." This thesis 1S intended to provide the 
designer a tool with which he can confidently accelerate 

the conceptual design phase and still produce high quality 
designs. The synthesis system will further free the designer 
to try innovative design concepts heretofor stymied by the 


burden of manual calculations and routine deciSions. 


ee SCOPE 

The purpose of this thesis is to couple two existing 
computer programs; a naval surface combatant synthesis 
model [{Ref. 3], and a general purpose non-linear optimizer 
[Ref. 4], to produce a synthesis system which will further 
enhance the ability of the naval architect during the 
conceptual design phase. Instead of using a synthesis model 
to generate hundreds of designs and then manually selecting 
one which appears to be the "best,'' the synthesis system 
will let the computer make the decision based on the limi- 
tations or constraints, and design requirements, coded in 
mathematical terms. 

This process may be illustrated using the traditional 
design spiral. Synthesis models currently in use reflect a 
computerized version of the spiral which iterate until they 
converge to a feasible design or diverge and terminate with 
Memaesion. Figure 1.1 illustrates this concept. Figure 1.2 
presents the optimizer/synthesis cycle proposed in this 


thesis. Here only one cycle is generated in the spiral with 


1s) 
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Figure 1.2 Optimized Design Spiral. 
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the optimizer changing design variables to meet the require- 
ments of the design simultaneously With the caiculations 
around the path. Differences between the desired and 
calculated values of specified parameters are treated as 
constraints which the optimizer works to satisfy. As can be 
seen graphically, much time and computational effort has 
been saved as a result of this process. 

It must be emphatically stated at this point that the 
synthesis system is a design tool, as are the individual 
synthesis model and the optimizer routine, and all of them 
must be used carefully with considerable common sense and 
good engineering judgement. 

The benefits of such a system are manifold and represent 
some of the best aspects of the individual programs. 

Several are listed below. 

1. An automated, optimized synthesis system further 
reduces the computational time necessary to do a 
meaoTOglauvy Study. 

2. The results are consistent throughout the design 
with the computer making decisions and calculations 
as a matter of routine. This consistency is impor- 
tant when considering designs of radical character. 

3. The results are more comprehensive. The synthesis 
model produces detailed output of essential data 
While the optimizer generates an easily traceable 
optimization trail, from which any design other than 
the optimum may be selected and evaluated. 

4. The ability to conduct studies on designs optimized 
with respect to different design objectives, i.e. 
minimum displacement versus minimum cost, while 


maintaining the same repeatable deSign practices 
and standards. 


A) 





5. And finally, the freedom to be innovative and 
creative in design as a result of no longer 
being tied to long, tedious manual calculations. 

D. PRESENTATION OF THESIS 

The remaining chapters of this thesis are outlined as 
follows. Chapter II presents a brief background on the 
optimization methods used in the COPES/CONMIN optimizer. 

A two variable design example is presented as an illus- 
tration to further highlight the optimization procedures. 

Chapter III is a summary of the Reed synthesis model's 
salient features including logical program flow and a 
brief description of the program's subroutines. 

Chapter IV presents the modifications that were necessary 
in the Reed model to couple it with the COPES/CONMIN opti- 
mizer. Constraints that were added are justified and 
presented. 

Chapter V is a short design example uSing mission 
requirements Similar to those of an FFG-7 as a baseline 
design. The example is used to present arguments sup- 
porting selection of objective functions and design 
variables. 

Chapter VI presents design examples using different 
objective functions and comments on the merits of the 
synthesis system. 


Chapter VII offers conclusions and recommendations. 


Oe 





Appendix A is a computer listing of the design example 
presented in Chapter V. 

Appendix B is a listing of the elements in the GLOBCM 
COMMON block, identifying their global locations and 


meanings. 


Appendix C is a listing of nomenclature used in this 


thesis. 
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II. OPTIMIZATION TECHNIQUES 


eee) LNT RODUCTION 

This chapter will present some of the fundamental 
concepts and definitions required to understand the opti- 
mization methods used in the COPES/CONMIN optimization 
program. COPES/CONMIN is then discussed in brief anda 
short two variable design example is presented to illustrate 
the methods used in determining the optimum solution. The 
explanations are necessarily short as it is assumed that the 
reader is already familiar with optimization methods. 
Should further details be desired, Fox [Ref. 5] and 
Himmelblau [{Ref. 6] are very good texts which may be 


consulted. 


me) DEFINITIONS. 

It 1S important to define several key terms used in the 
optimization problem formulation. These are: 

Design variables--The parameters for which values are to 
be chosen in vroducing a design will be called design vari- 
ables. In the conceptual ship design phase these might be 
length, speed, beam to draft ratio, length to beam ratio, 
block coefficient, etc. Design variables may be constrained 
to a limited range, i.e. 300 to 700 feet for the length 


oetween perpendiculars, or they may take on any value. 





Objective function--The single valued function with 
respect to which the design is optimized is called the objec- 
tive function. Selection of the objective function can be 
one of the most important decisions in the design process. 
For ship designs, varameters such as displacement, economic 
measures of merit, speed, structural member fabrication costs 
and weight may be used. It is also possible to combine 
several important parameters to form a weighted combination 
which is to be optimized. This procedure is discussed 
extensively by Leopold and Reuter [Ref. 7]. 

Constraints--The design restrictions which must be 
satisfied in order to produce an acceptable design are 
Sollectively called constraints. In the design of a ship, 
many constraints are specified by the owner prior to any 
calculations being done. Some of these may be maximum navi- 
gational draft, maximum or minimum length or displacement, 
and certainly some measure of stability. If a parameter is 
beyond the value of a specified constraint, the constraint 
is said to be violated. 

Infeasible Design--A design in which constraints are 
violated is called an infeasible or unacceptable design. 

Feasible Design--A design which satisfies the specified 
constraints is called a feasible or acceptable design. 

Side Constraint--A constraint which restricts the range 


of a design variable for reasons other than the direct 
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consideration of performance is called a Side constraint. A 
Side constraint in the design of a ship may be the minimum 


MWensth or draft. 


C. COPES/CONMIN 

The COPES/CONMIN optimization program is a general pur- 
pose, non-linear optimizer capable of handling large, con- 
strained problems. It has been successfully used in 
connection with aircraft synthesis models [Ref. 8], 
fervetural Optimization {Ref. 9], airfoil design [Ref. 10], 
and numerous other engineering applications. Although the 
code itself is very sophisticated, the primary methods of 
SaGimiZacion used are conjugate directions for locally 
unconstrained problems and feasible directions for locally 
constrained problems. 

It solves the general non-linear optimization problem 


stated as follows: 


Minimize OBS = CX) PN 
Subject to; G; (CX) < O iho <i = hm Eee 
ne ee For 2 = ln 2.3 

aL i a 


peere OBJ is the objective function. The vector X contains 
the "n" design variables. G,(X) define the constraints 
imposed by the designer on the optimization problem and "'m" 
is the total number of constraints. F(X) and G,(X) may be 
elther implicit or explicit functions of the design variables 


X, but must be continuous. The variables 5 and os Getine 
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the lower and upper bounds, respectively, on the design 
variable X, and are the limits over which F(X) and G, (X) 
are defined. 

The n-dimensional space spanned by the design variables 
Pees) terzerred to as the design space. As stated previously, 
any design which satisfies the inequalities of equation 2.2 
is referred to as a feasible design. If the design violates 
one or more of the inequalities, it is said to be infeasible. 
The minimum feasible design is said to be optimal. 

(ier Cwtatniazab ton program begins With an initial X vector 
which is input to the program and may or may not define a 
feasible design. The optimization process then proceeds 


iteratively by the following recursive relationship: 


(atl) _ x64) mn axg 64? 2.4 


mere Gg 145 the 1teration number, vector S is the direction 
of search in the n-dimensional design space, and a* is a 
scalar which defines the distance of travel in the direction 
Se 

The optimization process then proceeds in two steps. 
Miesiirst is the determination of a direction S$ which will 
reduce the objective function without violating constraints. 
The second is the determination of the scalar a* so that 
the objective function iS minimized in this direction, a 
new constraint is encountered, or a currently active con- 


traint is encountered again. 
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Consider for example, a hypothetical problem in which we 
wish to minimize the displacement of a naval combatant in 
terms of two design variables; length between perpendiculars 
(LBP), and the length to beam ratio (L/B). There are also 
two constraints for this problem: constraint Gl requires 
that the calculated full load disvlacement is equal to or 
greater than the displacement estimated by an empirical 
memmula, and constraint G2 requires that, in a similar 
manner, the calculated vertical center of gravity be 
greater than or equal to the empirical estimation. 

Figure 2.1 is a graphical representation of such a 
problem showing the contours of constant objective function 
value as well as the constraint boundaries. Assume that a 
@esign at point A is given so that initially no constraints 
are active or violated. The program then begins by 
ipectally perturbing each of the XK variables to determine 
its effect on the objective function (full load displace- 
ment). That is, the gradient of the displacement function 
1s calculated by the finite difference method. Because no 
constraints are violated or active, it is obvious that the 
greatest improvement in the objective function is obtained 
by moving in the negative gradient or steepest descent 
aeeeection SO that S = -V(DISPFL). Having determined 5S, 
the scalar a* in equation 2.4 must now be determined so 


that either the objective function is minimized in this 
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direction or some constraint boundary is encountered. A 
one-dimensional search is done in the direction $ to deter- 
mine the appropriate value for a* so that an improved design 
can be achieved at point B. No further improvement can be 
achieved in this direction as the constraint Gl is encoun- 
tered. Now a search direction must be found which will 
reduce the objective function without violating the 
@rsolacement constraint, Gl. Sweh a direction can be 

found by solving a linear programming subproblem with a 
Single quadratic constraint. The details of such a problem 
mee Fiven in [Ref. 9] and [Ref. 11]. 

The goal of the subproblem is to find a search direction 
which will minimize the objective and yet not violate any 
Semstraints. If no such direction exists, then the current 
point is considered optimal or at least a local minimum. In 
the example, a direction can be found and the design pro- 
ceeds from point B to point C where constraint G2 is 
encountered. The subprogram is solved again, resulting 
in a further reduction of the objective and an active 
constraint at point D. From point D, the one-dimensional 
Search yields a solution at point E in Figure 2.1 which 
is the vertex of the two constraints. Once again the 
design variables are perturbed to obtain the gradient of 
the objective and both active constraints, thus the linear 


subproblem is solved. This time the solution is zero, 





midtcatvane whe ODtimal deSign has been achieved. Point E is 
clearly the optimum since no re cision exists at this point 
which will reduce the objective function any further without 
violating one or both of the constraints. 

For preliminary designs, it is often quite possible for 
the initial design to start in the infeasible region. Logic 
meeinecluded in the optimization program so that if this 
BeeettaciOn occurs, a direction vector S is obtained which 
will point toward the feasible region with a minimal 
increase in the objective function. 

The methods used in this example are, in principle, 
directly extendable to the n-dimensional problem. Also, 
additional constraints can be imposed without increasing the 
complexity of the design process. For the Reed/COPES/CONMIN 
system it appears that the optimizer is comprehensive enough 
to easily handle all the design variables and constraints 
the designer could possibly want without any significant 
increase in computational time. Also, ship designs are 
merea for their condition of flat-laxity in terms of similar 
Ships having similar dimensions and design parameters. This 
SonmaLtion increases the probability of the optimizer 


achieving an optimal solution. 
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Pep eouee sor NIA eots MODELS 


ae HISTORY 

As stated in Chapter I, a synthesis model is an engi- 
neering procedure which converts a specific set of perform- 
ance requirements into the physical description of a ship 
which satisfies those requirements. There have been 
several synthesis models available for the naval ship 
designer since the late 1960's. Two of these models, 
currently used by the Navy, were the primary reference 
sources for the development of the Reed model used in this 
thesis. They are the U.S. Navy's destroyer model DDO7 
[Ref. 2] and Center for Naval Analysis Conceptual Design of 
Ships Model (CODESHIP) [Ref. 3]. A third model developed 
by the Coast Guard patterned primarily after the DDO7 
synthesis model was also used as a reference in constructing 
the Reed model [Ref. 12]. 

The destroyer model, DDOV, was first developed in the 
1960's and has been continually updated since then. This 
model was developed to represent only U.S. destroyer type 
Ships. 

The general framework for the method of analysis used in 
the CODESHIP model was also conceived in the sixties at the 
Center for Naval Analysis. CODESHIP was created to model 


Shaps which range in size from patrol craft to aircraft 


i 





carriers. Because of the wide range of ship types in the 
data base, the CODESHIP model is not as accurate in pre- 
dicting actual values as a model designed for a specific 
ship type, e.g. DDO? for surface combatants. The greatest 
asset CODESHIP provided in the development of the Reed model 


was its highly versatile input and output features. 


B. REED MODEL DESCRIPTION 

The Reed model used in this thesis was developed to pro- 
vide a more versatile synthesis model than CODESHIP and DDO’. 
Soecifically, the model allows the designer to control the 
design standards and practices to which the ship is designed 
ana to observe the resulting ship characteristics from a 
more functional level than existing models provided. Although 
based on standard calculations derived from functionally 
Similar ships, the designer is still given the opportunity 
to be flexible in his selection of values reflecting 
different design standards and practices. The resulting 
synthesis model enables the user to explore a significantly 
greater number of design options and provides a greater free- 
dom than previously allowed with manual calculations or even 
other synthesis models. 

The Reed model has been used at the Massachusetts Insti- 
tute of Technology since 1976. The primary purpose for the 
use there has been academic, although the results the 


model produces are certainly consistent with current design 
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peaeurces, In addition to its applications during concep- 
tual design, the Reed model can also be used during later 
design phases as a design aid for updating predictions, 

and as a tool for conducting comparative naval architecture 
studies. 

The Reed model is limited to surface displacement ships 
configured as naval combatants. The data base was created 
Means Ships which ranged in full load displacement from 1,770 
to 16,294 tons and which had lengths between perpendiculars 
of 301 to 700 feet. No cost characteristics are calculated 


by the model. 


C. REED MODEL ORGANIZATION 

The Reed model, as it was received, is controlled by a 
main program. The main program calls the various subroutines 
as required to calculate information for output and to let 
program execution proceed. Program control organization 
is shown in Figure 3.1. 

Mine logical program organization or program flow is 
shown in Figure 3.2. The program begins by reading the 
following data: 

1. Names of items to be printed in the output, 


2. Residual resistance coefficient values taken from 
the Taylor Standard Series [Ref. 13], and 


Eee Data for all the items in the payload shopping list. 
If more than one ship is to be run, the specifications for 


the second ship are then read using subroutine DATA2. 
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DATA2 - 


ca EPLANT 


MACHLQ MBSIZE 


Ea 
WEIGHT VRTCG 


SEASPD FNCGRP 
: OUTPUT 
VOLUME HPCALC 


SHEER CRVAL 


Figure 3.1 Reed Model Program Organization. 
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Following this, the main program calls subroutine SPAYLD 
which is used to sort the weights of the payload items input 
and place them in the proper Bsci? weight groups. The value 
of WPAYIN, the total weight of payload input, is calculated 
and is used for the initial estimation of the full load 
displacement, should the length between perpendiculars (LBP) 
not be specified. 

Next, the underwater hull shape is determined in the 
subroutine UWDIM. Values input to this subroutine are ep. - 
free surface effect correction, Cp, Cx, GM/B, as well as the 
estimated KG and full load displacement. Output values from 
UWDIM are beam, draft, and Cwp. 

If the length between perpendiculars LBP, is not speci- 
meed at the start of the program, it is calculated in the 
subroutine GEOM prior to transferring control to UWDIM. 

This feature was removed from the synthesis model when it 
was coupled to COPES/CONMIN and it is assumed that the 
designer will have a reasonable estimate of LBP to begin 
the design. 

Subroutine HPCALC is called to determine either the 
maximum sustained speed, VSUS, or the horsepower required at 


13SCI--"Bureau of Ships Consolidated Index of Drawings, 


Materials, and Services, Related to Construction and Con- 
version" of February 1965 and given in an appendix to 
Reference [14]. 


“variables used in this thesis are described in Appendix 
C, Nomenclature. 
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the maximum sustained speed, SHP. If VSUS is input, the 
HPCALC routine along with the CRVAL subroutine generate an 
estimate of the horsepower required to maintain the maximum 
Sustained speed. If SHP is input, HPCALC is used to compute 
a corresponding VSUS through an iterative procedure. HPCALC 
is called a second time to find the horsepower required at 
endurance speed. Endurance speed must be input. 

Next, the subroutine EPLANT is called to calculate 
cruise, battle, and 24 hour average electric loads. Based 
on these estimates, the required number and size of the 
generators is determined. The electric loads and number of 
generators may be input, in which case these calculations 
will be by-passed. 

Subroutine MACHLQ is called to calculate machinery and 
other related liquid weights. Liquids considered machinery 
related are potable water, reserve feed water, lube oil, 
fuel oil, and diesel oil. These calculations are based on 
the military requirements of the ship. Liquid weights are 
used to determine the required tankage volume. 

MBSIZE determines the minimum depth of the machinery box 
which then corresponds to the minimum depth the ship must 
have amidships. 

The next step is to calculate the required and available 
volumes for the ship and to achieve a balance between the 


two. This is done in the VOLUME and SHEER subroutines. As 
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most naval combatants are now considered volume-limited, 
this is one of the most important aspects of the synthesis 
model. 

The subroutine WEIGHT next determines the weights of all 
light ship and load items that are not payload-related. 
Full load displacement is calculated and it is then compared 
to the estimated value of displacement used to proceed this 
far in the program. If the two values are not within a 
specified tolerance, say 10 tons, the program reestimates 
displacement and iterates through all previous calculations. 
If the two values are within the specified tolerance, then 
the program continues on and compares the estimated and 
calculated values for the vertical center of gravity. Ina 
Similar manner, if the two values are within a specified 
tolerance, the program continues; if not, there is a new 
estimate generated and an iterative process proceeds until 
agreement is made. The vertical center of gravity is 
calculated in the subroutine VRTCG. If after a specified 
number of iterations, the comparisons for either displace- 
ment or vertical center of gravity are still not within the 
desired tolerances, the program prints an error message and 
terminates. 

At this point in the synthesis model, all essential 
calculations have been completed and the principal charac- 
teristics and coefficients of the "feasible" design have 


been determined. FNCGRP is called to place weights and 
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volumes into their proper functional groupings as listed in 
Appendix C of Reference [15]. These groupings form a major 
portion of the detailed output from the synthesis mcdel. 

The subroutine SEASPD calculates the speed the vessel is 
able to sustain in the North Atlantic Ocean in order that 
relative comparisons may be made of seakeeping characteris- 
mcS . 

imewitinal subroutine 1s OUTPUT which prints the input 
specifications first, followed by the payload items input, 
summary of results from the final estimate of principal 
characteristics, functional group results, BSCI weight 
listings, and the functional electrical loads. The OUTPUT 
routine is very versatile in that any portion of the above 
output may be supressed or printed as the designer desires. 

The above descriptions of the subroutines in the Reed 
Synthesis model were designed to give the reader a general 
idea of how the program functions. The routines are 
Meccribed in much greater detail in [Ref. 3] and [Ref. 12]. 
The reader is encouraged to use these to make changes to the 
synthesis model or to gain a greater understanding of the 


program. 
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IV. COUPLING COPES/CONMIN AND THE REED MODEL 


A. SELECTION OF DESIGN VARIABLES AND OBJECTIVE FUNCTIONS 
One of the most important decisions in the development 
of the synthesis system is the choice of parameters to be 
used as design variables and objective functions. Review of 
current literature pertaining to preliminary design syn- 
thesis and optimization models, [Ref. 16-22], produced two 
different approaches to the problem. Nowacki [Ref. 16], 
worked with a single economic measure of merit (required 
freight rate) for the objective function of tanker designs 
and used speed/length ratio, B/T, L/B, Cb, and L/D as 
design variables. Mandel and Leopold [Ref. 17], recommended 
that for tanker or general cargo vessels a three-term 
weighted optimization criterion be used as the objective 
function. One term would be an economic criterion, the 
second term would take into account payload weight and the 
third term payload volume. Both payload terms are based on 
the owner's requirements. They selected displacement, Cp, 
speed/length ratio, B/T and L/D as the five design variables. 
Leopold and Reuter [Ref. 7], carried the multiple term 
Optimization criterion a step further in their paper on 
design methods and philosophies for the FDL, LHA, and DD-963. 
Here they proposed an optimization criterion containing 


seven terms; cost, effectiveness, flexibility, availability, 
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naomi ptlaey,) Vulnerability, Survivability, innovation, and 
nonobsolesence. Through this was never expressed as a 
mathematical function, it indicates the attempts to 
quantify key parameters in optimization methods. 

From the summary presented by Nowacki [Ref. 16], it was 
apparent that some sort of economic measure of merit was the 
dominant objective function for commercial vessels. This is 
to be expected as profit is the motivation for construction 
and operation of a commercial vessel. Military vessels, 
specifically naval combatants, are not required to be profit- 
able. Their mission is to deliver a military payload at 
a specified time and place and to provide military services 
when and where needed. Therefore, economic criteria, though 
certainly important, are not as vital to the designer as are 
the ability of the vessel to carry its military payload and 
fa—rorm its military duties. 

In recent years, the naval combatant has become area/vol- 
ume critical in nature. Accordingly, this has become the 
dominant factor in establishing the ship size. Taking into 
account the fact that volume is critical and that there 
exists an implicit relationshivd between volume and displace- 
ment, it was decided that displacement would be the most 
representative variable for the objective function. Accord- 
ing to Manning [Ref. 22: p. 101], using least displacement 
aS a measure of merit for the military vessel gives the 


highest ratio of military payload to displacement. This 


of 





selection is further supported by the fact that total ship 
cost can be roughly estimated by the displacement. 

The design variables selected for the Reed/COPES/CONMIN 
system were LBP, L/B, B/H, Cp, and Cx. These variables were 
selected by running sensitivity studies on displacement with 
respect to each design variable over a specified range of 
values. The payload and all other factors remained constant 
in this analysis as each design variable was tested. It is 
felt that these variables are the most significant in the 
Gerermination of the principal characteristics and provide 
the designer with considerable flexibility for innovation. 

It should be noted that any one of the design variables 
may be designated as the objective function so that compari- 
sons of designs with, say displacement as the objective 
function can be made with designs that may have used LBP or 
Cp as the objective function. Also, should the designer 
desire, he may designate any of the other variables used in 
the Reed model as objective functions and/or design vari- 
ables. It is only necessary to ensure that the variable is 


listed in the GLOBCM statement. 


B. THE GLOBCM STATEMENT 

The GLOBCM COMMON statement is a requirement of the 
COPES/CONMIN program. All variables used as objective func- 
tions, constraints and/or design variables must be listed 


in the common statement and the statement must appear in 
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each subroutine the variable is used in. It is used by the 
optimizer as a catalog to identify where the design vari- 
ables, objective function(s), and constraints are, and what 
purpose they fulfill. In the Reed/COPES/CONMIN system the 
GLOBCM statement has DISPFL, full load displacement, and a 
meme leonent vector labelled S as entries. The 5S vector 
contains input parameters, payload items and quantities, 
BSCI weight group data, BSCI vertical center of gravity 
data, volume values for payload and input items, elec- 
trical load values, miscellaneous parameters and coefficients, 
and space for inputting any special payload items not in- 
cluded in the standard payload shopping list. DISPFL was 
meed tO GLOBCM because it was not in the § vector. All of 
the design variables used in this thesis are contained in 


fees vector. The element numbers are listed in Appendix B. 


C. CONSTRAINTS 

In order to couple the Reed model with COPES/CONMIN, it 
was necessary to make several major modifications to the 
synthesis program. [In addition to the modifications, all 
default parameters were replaced by constraints. 

The main program in the Reed model had two iterative 
loops which were removed. The first loop occurred after all 
the weight calculations had been done. The purpose of this 
loop was to check to see if the calculated value for the full 


load displacement equalled the estimated value for displacement 


= he, 





within a user specified tolerance. [f the difference was not 
within the tolerance, the program would estimate a new 
value for displacement, calculate a new full load displace- 
ment and make the check eon The iteration continued until 
the check was satisfied or a maximum number of iterations 
was achieved at which time an error message was printed and 
the program terminated. This loop was replaced with a single 
constraint. The quotient of the calculated displacement to 
the estimated displacement minus one is required to be less 
uoam Or equal to zero. 
G1 = DISPFL/DPTRY - 1.0 < 0.0 

The second loop in the main program is exactly the same 
as the first except that it compares the initial and cal- 
culated values for the vertical center of gravity (VCG). 
This loop was replaced with a constraint similar to the 
first using calculated VCG and estimated VCG. 

G25 =" CGriD/ KGPRY = 1.0 < 0.0 

In the subroutine UWDIM, two default values were replaced 
by constraints. These defaults required that the beam to 
draft ratio be greater than two and less than four in order 
to use the residual resistance coefficients in the Taylor 


Standard Series tables. The constraints are: 


Gem= 2.0/8TR = 1.0 < 0.0 


G4 Bon; 4.0 — 1.0 5s 0.0 


i 


A third constraint is added here to satisfy a stability 
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Pere rleote Aevarraole HR 1S introduced in the subroutine 
with a value of: 

R = KB + BM - GM - KG - FSCORR 
When R = O, the stability criteria is just satisfied, there- 
fore the constraint 

GS = -R/10 < 0.0 
must be less than or equal to zero to ensure a stable de- 
Sign. R is divided by ten as a means of scaling in order 
that GO will have the same relative Seer as the other 
constraints. 

Subrovtines HPCALC and CRVAL have the greatest number of 
constraints in them due to the limitations of the Taylor 
Standard Series power estimation. The maximum ranges for 
the input variables that are accepted are listed below: 

ye) Sey ce. 0 
OMe CO <) 0.70 


WeOOL = Cy = 0.006 for 0.5 


1A 


iy Iie feces 


oot =< Cy =< 0.005 for O.5 


1A 


WO 56, Se 2 a8 


These are translated into the following constraints: 


G3 = 2.0/BTR - 1.0 < 0.0 
G4 = Bley 4.0 = 1.0 < 0.0 
Cam atop — £.0 < O40 
Gq = €p/O.70 = 1.0 < 0.0 
Goe=w0eo/ SURAT = 1.0 < 0,0 
Meer moter /1.3 - 1.0 < 0.0 
ecu OnOOny ey = 1.0 < 0.0 
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GiG—sev, 8,006 — 1.0 < 0.0 

Gilli = Sevat/2.0— 1-0 < 0.0 

POF Cy/ OmbOs - Laoins O50 

G12 = -1.0 

iG GenclmeOnMOD mand. (XMis.Gl. OF0)) ~GilZ = 1.0 

The remaining constraints appear in the subroutines 

VOLUME and SHEER. In the routine VOLUME, there was a loop 
which iterated to match the total estimated ship volume with 
the total calculated ship volume by increasing the deck 
house size or by adding a raised deck to the ship. This 
loop was replaced by the following constraint: 


ioe —mnoov/DEHY¥e= 1-0 <.0.0 


D. FURTHER PROGRAM MODIFICATIONS 

As stated previously, it is necessary for the operation 
of COPES/CONMIN to have the analysis portion (i.e. the 
synthesis model) of the program in subroutine form. The 
Subroutine is called ANALIZ and has the calling parameter 
ICALC. When ICALC = 1, the data for the operation of 
SerPrs/CONMIN is read in. For ICALC = 2, all the analysis 
calculations are performed by the optimizer, and when [CALC 
= 3, the final results are printed out. Lines were inserted 
into the Reed program at the appropriate points to facilitate 
the operation of ANALIZ as stated above. 

In the stand alone form, the Reed model is capable of 


analyzing one or more ships as the designer desires. The 
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Reed/COPES/CONMIN system will process only one design at a 
time. This is not considered a problem because changes-to 
the input data files are made easily and the computational 
meme fOr each run is short. 

Subroutine GEOM was deleted from the program. This 
subroutine calculated a LBP based on the weight of the 
payload when no LBP was input to the program. It is felt 
that the designer can estimate a LBP which the optimizer 
can start with and proceed to an optimum design. 

Extensive changes were made in the subroutine SHEER. 

It is assumed that the minimum depth amidships is the same 
or greater than the depth of the machinery box. This 
value is compared to LBP/16 and the largest of the two 
values is selected as the depth amidships. Having satisfied 
this criterion, an iterative routine was used to estimate 
the freeboard at the forward and after perpendiculars. 
The iterations were based on an assumption that sheer 
Meaecctions Of 0.01 and 0.03 forward and 0.001 and 0.0075 aft 
were esthetically and structurally satisfactory. If these 
fractions could not be met, the program printed an error 
message and terminated. This was changed so that when the 
depth amidships was established, the freeboard at the forward 
and after perpendiculars were assigned as: 

ZO wCOnOStiBE + DLO) = HA 


F20 = (0.0075*LBP + D100) - H 
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respectively. These values are comparable to those produced 
by the Reed model alone and simplify the design process. 
Should different values be desired, it is only necessary to 
change the two equations. 

The Reed model was written giving the program the option 
moenadd a raised deck to the hull if the total required 
arrangements volume was greater than the total available 
arrangements volume. The option was deleted from the 
program and replaced by constraint G13. Should the designer 
wish to have a raised deck design, the calculations can be 
replaced but care must be taken to ensure that no design 
variables are altered by an internal loop, thereby short- 
cutting the optimizer. 

There were numerous small changes made to the Reed model 
to make the coupling of the two programs more efficient. 
They are, however, insignificant and will not be documented 


here. 
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V. DESIGN EXAMPLE 


fe INTRODUCTION 

This chapter will present an example using the 
Reed/COPES/CONMIN synthesis system. It is assumed that the 
user is familiar with Appendix A, the User's Manual, of 
Reference [3] and the data description and input format for 
COPES presented in Reference [4]. The logic for the selec- 
tion of the objective function, design variables, and 


constraints will be discussed as the method of problem 


presentation. 


B. EXAMPLE DESCRIPTION 

All designs start with the owner's requirements or in 
military terms, the Operational Requirements (OR). This 
delineates the military payload and mission to be accom- 
plished over the lifetime of the ship. From the OR, the 
designer is able to select those parameters which will 
help him make an initial estimate of the ship, i.e. start 
the conceptual design phase. In this example, the 
operational requirements will be similar to those of an 
AAW/ASW capable escort ship like the Oliver Hazard Perry 
class frigate. 

Mission requirements specified in the OR will generally 
dictate the type of equipment required in the design. Table 


21 of Appendix A, [{Ref. 3] can be used like a shopping list 
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at this point. The designer can pick and choose the weanons 
systems, electronics suites, ammunition load-outs, prelimi- 
nary liquid loads, and aviation items he desires to use. 
This list of items technically defines the military payload 
hardware. Additional men and materials are required to 
support: and operate the payload. A partial list of payload 
items selected for this design example are listed in Table 
if. 
TABLE I 
Military Payload List 

FFCASW command and control) Radio communications 

SPS-49 Radar w/IFF 

SPS-55 Radar w/IFF 

SQS-56 Sonar 

FF/FFG Basic ECM suite 

ASWC and C-FF-2C,7D Electronic Tactical Data Systems 

Vulcan/Phalanx on O1Lv 

76mm Gun, OlLv 

Mk-92 CIWS/STIR 

evO73''/90' rounds 

10000 20mm rounds 

Mk-13 Tartar Missile Launcher w/40 Missiles 

Mk-32 Triple Torpedo Tubes P/S w/24 Torpedos 

Harpoon FCS 


1 Lamos MK III Helo with support 
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Having made these decisions, which incidentally may have 
taken many, many weeks, the designer then proceeds to select 
those parameters which define the ship form. Here the 
designer has slightly more freedom than with the payload 
entries. It is possible to run systems trade-off studies 
comparing different propulsion plants, hull materials, and 
design standards. Items in this category come from Table 12 
of Appendix A, Reference [3]. They are such things as type 
of propulsion plant desired, number of propellers, the 
endurance range and speed, geometric parameters such as 
Memeth, L/B, B/H, coefficients of form, electric plant 
parameters, and crew accommodations. 

Up to this voint in the process the synthesis model and 
the Reed/COPES/CONMIN system are essentially the same. To 
continue further with the synthesis model will require that 
relatively accurate values of the ship specifications be 
input. The Reed/COPES/CONMIN system on the other hand 
allows any starting value that is within the specified 
limits of the design variables. That is, it is possible to 
start with an entirely infeasible design and yet end up with 
a feasible design after optimization. This is one of the 
strengths of the synthesis system. 

From the specifications identified in this group of 
parameters, the designer selects the objective function(s). 


design variables, and constraints. He may also develop a 
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function he considers viable as an objective function and 


insert it into the synthesis model. Table II lists those 


variables which must be invut in the problem specifications. 


TABLE II 


Shiv Specifications List 


Sustained Speed 
Endurance Speed 
Endurance Range 

LBP 

L/B 

By Hi 

Cp 

Gx 

Propulsion Pit Type 
Sustained SHP 
Number of Main Engines 
Number of Shafts 
Propeller Type 


Shaft Type 


Ship Service Elect Plt Type 


C. PARAMETER DESCRIPTIONS 


Emergency Electrical Plant Type 
Type of ship heating 

Fin Stabilizers (Y/N) 
Officers' Accommodations 

CPO Accommodations 

Enlisted Accommodations 

Flag Accommodations 

Troop Accommodations 
Passenger Accommodations 
Basic Hull Material 

Basic Superstructure Material 
GM/B Stability Value 

Design Margin 

Free Surface Correction 


Passageway Type (P/S or CL) 


Selection of the objective function in this example is 


based on the comments of Manning in Reference [22] and 


Saunders in Reference [23]. 


From Manning: 


"There are two criteria which may be used to measure 


the excellence of a design. 


The one usually applied 
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to warships is that the best design is the one of 

least displacement. This design gives the highest 

Vullemoet themiatio on military load to displacement." 
Accordingly, it was decided to use displacement (DISPFL) as 
the objective function. The intent will be to minimize the 
displacement required to support the specified payload items 
and to use this as an indicator of relative size and cost 
when compared with other designs. DISPFL is listed in 
global location one in the GLOBCM statement and identified 
as the objective function to be minimized in data block E. 

Five parameters have been selected as independent design 
variables. The first is the length between perpendiculars. 
Justification for this selection is well stated by Saunders 
on page 343 of Reference [23]. 

erm the group of underwater form coefficients and 

parameters developed through the years, the ship 

length L logically appears as one of the principal 

dimensions. It is related directly and indirectly 

to the beam, the draft, the displacement weight, 

the displacement volume, and to many other factors." 
Length is designated by LBP and is in global location 5. 

The second design variable is the ratio of length to 
beam (L/B). Although the L/B ratio can be related to the 
prismatic coefficient, Cp, and the B/H ratio and is there- 
fore not strictly independent, it has still been selected as 
a design variable. Motivation for this selection comes from 
the methods used in the Reed synthesis model to determine the 


beam and initial powering estimates. It may also be used as 


a weak measure of turning characteristics [Ref. 23: p. 352] 
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and course tracking when comparing designs. The L/B ratio 
imme clobal location 6. 

The beam to draft ratio (B/H) is the third design vari- 
able. This ratio is used to determine the ship's draft and 
is used in estimating horsepower requirements. It is stored 
mame tobal location 7. 

The fourth design variable is the prismatic coefficient, 
Cp. This coefficient indicates the fullness of the under- 
water hull and is therefore indirectly related to displace- 
ment. A small Cp means the ship has fine ends anda 
full midbody. A large Cp means full or blunt ends similar 
moma canker or barge. Cp is in global location &. 

The last design variable is the midship section coeffi- 
cient. This relates the area of the midship section to a 
rectangle whose sides are equal to the draft and the beam at 
that section. It is used in estimating the hull strength 
and initial power requirements. It is located in global 
fmecation 9. 

With the design variables identified it is now necessary 
to specify their global location in the GLOBCM statement and 
to specify any side constraints that may be imposed on the 
design variables. COPES data block F is used to specify 
Side constraints and data block G provides COPES with the 
design variable number and global location. Design vari- 


ables i-5 are located in global locations 5-9. Table III 





presents the side constraints that are imposed as a result 
of the restricted data base used in the synthesis system. 
TABLE III 


Side Constraints on Design Variables 


PARAMETER LOWER BOUND UPPER BOUND 
LBP SOOmrt . OOM S - 

L/B 700 120 

B/H 2.00 4.00 

Cp 0.50 0.90 

Cx 0.75 0.90 


Bee CONSTRAINT PHILOSOPHY 

Meets Worch, co philosophize a bit at this point about 
side constraints and about constraints in general. The 
constraints imposed on the design variables in this example 
are used to maintain the design within the specifications of 
the synthesis model data base. Extrapolation of designs 
based on data and on empirical relationships derived for the 
model is considered risky at best and would probably be 
better accomplished with manual calculations. Further, the 
design is constrained by the limitations of the residual 
resistance coefficient values of the Taylor Standard Series 
used in the horsepower estimations. The obvious point to be 
made is that as the synthesis system is constrained, so the 
designer is constrained. He cannot develop feasibility 


Scuaqies for aircraft carriers or air cushion vehicles with a 
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destroyer synthesis model. Also, he cannot justifiably try 
radically different design concepts for a destroyer if the 
concepts render the data base invalid. The manner in which 
the designer attacks this problem is an indication of his 
creativity and foresight. Improvement in the synthesis 
system is certainly an option to help eliminate some of the 
constraints required. For this synthesis model, the resis- 
tance coefficients could be expanded through the testing of 
ship models whose size is larger than those used to obtain 
the Taylor Standard Series data. This is an expensive and 
time consuming process and must be considered carefully in 
light of what is desired from the synthesis system. Again, 
the designer's ability to use the synthesis system as a tool 
is the key to his success. 

Hypothetically, a synthesis system should be designed to 
model real world conditions and require only those constraints 
imposed by the real world, i.e. no negative areas or plate 
thicknesses. The Reed model is considered an excellent 
synthesis model and requires relatively few constraints to 
maintain real world conditions. In some problems, it is 
necessary to use constraints to prevent the optimizer from 
exploiting a weakness or undesirable trend in the analysis 
portion of the code. This is not the case with this system 
and this point is mentioned only in passing should the designer 


desire to make significant alterations to the synthesis model. 
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Howumewp rte lyveatscussed constraints, it 1S now appro- 
Dreiate to identify them specifically and illustrate how they 
are delineated in the COPES data cards. In Chapter IV 
thirteen constraints (G1-G13) were proposed to replace 
several iterative loops and termination criteria built into 
the original Reed model. All of these constraints are 
Satisfied when they have a value less than or equal to zero. 
These constraints were assigned unused numbers in the $ 
feemor and occupy global locations 77 through 89. All of 
the constraints together form a constraint set as they are 
all of the same magnitude and are satisfied by meeting the 
Same criteria. Data block H is used to specify the number 
ppemconustraint sets, in this case one. Data block I identi- 
fies the first and last global locations of the constraint 
set, 77 and 89, and specifies the upper and lower bounds and 
any scaling desired for the constraints. For this example, 
the upper bound is zero, the lower bound is numerically 
minus infinity and the scaling value is defaulted to 0.1. 

Should any further constraints be desired the designer 
weoowmouly 1dentity them in the S vector, specify their 
eiieeal locations, their upper and lower limits and any 
scaling necessary via the appropriate data blocks. Table 
mets a listimg of the COPES/CONMIN data blocks illustrating 


the entries described above. 
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VUES, OVE 
COPES/CONMIN Data Blocks 


$ BLOCK A 

OPTIMIZED VERSION OF THE FFG-7 
$ BLOCK B 

Ds 

ST BLOCK C 

3,,,10,,,90 

$ BLOCK D 

O75, 0-02 


© 


LBP 

SCO OMe 76.010 @. 

os 
a0 aah) 

B/H 
204.0 

Cp 
ro0-0.90 

Cx 
wo .0590 
BLOCK G 


PAU WNHORrRrMOWMOWMANYW+~)1%A 


Sebo. © sO nO ~O 0 
$ BLOCK V 
END 
E. EXAMPLE RESULTS 
The results of this example are best presented in 
tabular form. Table V illustrates the initial and the 


optimum values for the objective function, design variables, 


and constraints. Note that the initial design was infeasible 
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Waeieconceragints Gl, G2, and GlS3 violated. Physically, this 
meant that the initial design was unstable (G2) as a result 
of insufficient volume (G13) and displacement (G1) to support 
the required payload. The optimizer increased the displace- 
ment by manipulating the design variables until all the 
constraints were satisfied. Constraint G13 is indicated as 
active, which means that it lies within a specified tolerance 
value of the constraint zero value boundary. 

TABLE V 


Initial and Optimum Parameter Values 


PARAMETER INTERAL VALUE OPTIMUM VALUE 
Objective Function 
DISPFL 2865.05 3511.68 
Design Variables 
LBP 200:. 0 394.38 
L/B 9.07 8.76 
B/H 3.14 cee 
Cp 0.59 0.50° 
Cx Oras: OA 767) 
Senstraints 
G1 6222 - -0.015 
G2 eo" -0.068 
G3 Sos -0.398 
G4 Se -0.168 
G5 Hy) falls -2.062 
G6 =o @ -0.040 
G7 S52 -0.285 
Gs -~5.66 -0.503 
G9 35). Sy ~0.501 
G10 iS ON -0.665 
Gll BP) -0.496 
G12 ORO) -1.000 
G13 ie 86e -0.0006° 


‘lower bound of design variable 
*violated constraint 
Sactive constraint 





Figure 5.1 is a graph of the objective function versus the 
number of iterations required to reach the optimum. It 

can be seen that the optimizer is indeed efficient and in 
this example arrived at a value very near the optimum in 
five iterations. The remaining four iterations were done 
Simply to "fine tune" the value of the objective function. 
The average run time for this example was approximately two 
and one half CPU second on an IBM 3033 system. The entire 
example with values for all the parameters at each iteration 
is provided in Appendix A. 

Comparison of the data arrived at in this example using 
the general specifications of the FFG-7 were very good. 
Displacement was within five percent of the actual ship dis- 
placement. Most of the other parameters were within at least 
ten percent of the actual and could have been made more 


accurate by greater detail in the problem specifications. 
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Figure 5.1 Objective Function vs Iterations. 
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VI. SYNTHESIS SYSTEM APPLICATION 


A. DESIGN EXAMPLES USING DIFFERENT OBJECTIVE FUNCTIONS 

The Reed/COPES/CONMIN system is extremely versatile in 
application. The designer is able to optimize a design for 
a svecific payload with respect to a larger number of impor- 
tant and different parameters. One example has been shown 
using DISPFL as the objective. LBP could be used as an 
objective function. Based on the relationship between 
length and displacement, the designer would minimize length. 
Another parameter that could be used is SHP/TON. The 
designer would maximize this function in an effort to get 
the most horsepower in the smallest ship. Table VI presents 
several examples which were optimized with respect to 
different objective functions. The baseline ship is a close 
approximation to the FFG-7 ocean escort. The same design 
variables were used in all cases. The constraints were not 
altered and the payload was also maintained constant. Data 
produced by the original Reed synthesis model is included 
mer comparison. 

Several general comments can be made from an overall 
view of the data. First, is that there appears to be a 
rather well defined optimum ship regardless of the objective 
Mmemetlon. This Supports the so called "flatlaxity" of ship 


designs [Ref. 16: p. 105] which allows the optimum ship to 
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have slightly different principal dimensions yet remain 
optimum or very nearly so. The second comment is on the 
consistency of the synthesis system. When compared with the 
original synthesis model results, the optimized results 
accurately reflect the effect of changing the design vari- 
ables to satisfy the design requirements. Additionally, 
specific parameters demonstrate excellent repeatability from 
design to design. This indicates that the final design is 
mie ODtimum and not a relative maxXimum or minimum point in 
the design hyperspace. 

The above comments are well supported by the data in 
Table VII. This data represents four deSigns: one starting 
with LBP equal to 300 feet, one with LBP equal to 500 feet, 
one with LBP equal to 700 feet, and a fourth design using 
parameters from the optimum design achieved with LBP equal 
to 300 feet. All cases were run with DISPFL as the objec- 
tive function and all other parameters held constant run to 
run. Excellent repeatability and a well defined optimum are 
demonstrated. 

A third and interesting point to note is that not all 
selected objective functions or variables are affected by 
the ovutimization process. Note the value of VPAY/VOL 
(Volume Payload/Total Volume). To bring about a change in 
this variable, the synthesis model msut be changed. This 


represents a change in the design standards or the designer's 
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Bielosephy., “lhe ability to be able to make such alterations 
famine ind!cation Of the strength and flexibility of the 
Synthesis model. A particularly good example of another 
such change is with regard to habitability requirements, 
i.e. ft*’/man. Current designs have shown an increasingly 
large portion of internal volume devoted to habitability. 
It is possible through the use of scale factors to modify 
current design standards established for habitability 
and to determine the effect on the overall ship design of 
changes in these standards. An excellent example of this 
feature is presented in Reference [3] where an FF 1052 
class frigate is redeSigned using Soviet habitability 
standards and then compared to the actual frigate design. 
From Table VI, the design which deviates most from 
the optimum is the one in which the objective function was 
full load density (FLD DENS). Here the objective function 
was to be maximized. The units of full load density are 
lbs/ft’? and are therefore a measure of how dense the ship 
is. The difference between this design and the others is 
Significant and leads to a consideration of which is more 
eertelent from a space utilization standpoint. It would be 
a logical assumption that the larger ship could support a 
bigger payload more efficiently. Yet, is this the better 
design? The designer must remember that the synthesis system 
is a design tool and he must be judicious in the application 


of engineering judgement and common sense. Optimization 
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cannot be blindly applied to a problem and the results 
accepted at face value. Clearly, other considerations would 
be necessary to choose between this and one or more of the 
other designs. 

As with the example in Chapter V, all the designs began 
in the infeasible region and proceeded to a feasible design. 
This feature is of great benefit to the designer as numerous 
designs that the synthesis model would have rejected and 
required new information for are simply tried and identified 
as infeasible on the path to a feasible design. The syn- 
thesis system maintains a record of all the designs tried 
which the designer may scrutinize and use for future design 
decisions or historical data. The system also identifies 
those constraints and design variables that are active or 
violated, thus providing the designer information on what is 


Mesey Critical in the design. 


Eee SENSITIVITY STUDIES 

It is possible to generate sensitivity studies with the 
Reed/COPES/CONMIN svstem. Once the designer has identified 
a design he wants to work with, he is capable of determining 
the effect of varying a single parameter at a time while 
Maintaining all other parameters constant. This was one of 
the first forms of preliminary design optimization used and 
is described in the landmark paper by Murphy, Sabat, and 


Taylor in 1965 [Ref. 15]. This feature was especially 
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helpful in determining the design variables and objective 
functions to be used in this thesis. Figure 6.1 illustrates 
a simple example where the optimum design minimized with 
respect to displacement is checked for sensitivity to 
variation in length. Note that while displacement can be 
reduced further by making the LBP less than 394 feet, this 


Violates one or more constraints. 


C. SYSTEM TRADE-OFF STUDY 

An important and useful application of the synthesis 
system is systems level trade-off studies. The Reed/COPES/ 
CONMIN system allows these studies to be conducted with 
consistency and confidence as the computer makes routine 
decisions and calculations the same in all cases. Table 
VIII presents data from a systems level trade-off study 
done on different types of propulsion plants. The second 
ceneration gas turbine (2 LM2500's) appears to be the best 
design while the steam and diesel plants are all comparable. 
Surprisingly, the COGAS plant appears to be the worst of 
the designs studied. This may be because the data used 
to synthesize the COGAS design may be dated and needs to 
be revised. As in the other examples presented in the 
thesis, the payload requirements for these ships are Similar 
to an FFG-7 class frigate. 

There are considerably more applications available with 


the Reed/COPES/COMIN system which have not been explored. 
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Figure 6.1 Sensitivity Displacement to Length. 
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TABLE VIII 


Systems Trade-off Study on Propulsion Plants 


PARAMETER STEAM STEAM GAS DIESEL COGAS 
600 lb 1200 lb TURBINE PLANT PLANT 

De 417.49 418.86 394.33 416.84 468.74 
BEAM 03.81 Dose 46.89 47.70 09.64 
DRAFT 14.951 14.56 113) 5 24 13.47 Nlrs peel Ere: 
Cp One a 0.52 0.50 Or 0 Orsi 
Cx 0.80 0.80 O12 O90 OF 10 
VCG FLD IS) 5 Sie: jG) Elere! 13.40 1S 5 2b 15.20 
L/B tha TAS Tels 8.42 8.74 (ele: 
B/H 657k SiG (fl 3.46 3.54 8) Ae: 
Vsus 28.46 28.52 Bee) SP Be) Veo) ae oO 
DISPFL 4129 4164 S15) ial 4110 0958 

DISPLS 2991 3029 2471 3091 3294 

FLD DENS 1S) 5 fox, Gao eye al oo Zl G 
WPAY/FLD On09 O09 ORAL One? O20 7 
WPER/FLD O03 O10 3 0.04 O03 0.03 
VOLPAY /VOL Ono 0219 OG 28, 0.18 CEs 
VOLPERS /VOL O20 og AG O@22 Ong 0.20 
VMB/SHP Soe 3) 9 By 2.40 4.26 2.40 
SHP/DISPFL 9269 So el ioe 2. 4o eae 
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COPES/CONMIN is capable of performing optimum sensitivity 
Studies, two variable function Space studies, and approxi- 
mate optimization techniques. Additionally, numerous 
combinations of variables and design standards are available 
in the Reed model. The examples presented here and in the 
previous chapters just begin to illustrate the utility of 


such a design tool. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


The Reed/COPES/CONMIN system provides the naval archi- 
tect with an efficient design tool, with which he is capable 
of reducing the time required to perform feasibility and 
concept design studies. The designer has at his fingertips 
all the benefits of a sophisticated synthesis model which 
has been proven to produce reliable designs, as well as the 
power of a nonlinear numerical optimizer. Combined together, 
the synthesis/optimizer system enables the designer to be 
more exhaustive in his search for a design solution at no 
greater cost in time. A greater number of design alterna- 
tives can be vrocessed and compared as a result of the 
synthesis system's ability to start with an infeasible 
design and develop it into a feasible design. The designer 
is made aware of the design variables and constraints which 
are critical in each case and has the ability to alter those 
and other variables to achieve his goals. 

In addition to being computationally fast, the synthesis 
system is very versatile. The COPES/CONMIN optimizer is 
capable of doing optimization analysis, sensitivity studies, 
optimum sensitivity studies, and optimization using approxi- 
mation techniques. The Reed synthesis model is capable of 
accurate conceptual design calculations of greater detail 


than previously attempted in synthesis models. It allows 
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for changes in design standards and philosophies as well as 
equipments and characteristics. Studies conducted using 
various combinations of features available from both 
programs greatly enhance the information available to the 
designer at any time. 

There are several areas where further development would 
be desirable. The first is the addition of a cost estimating 
routine. This module could include lead ship and follow 
on ship acquisition costs as well as life cycle cost 
analysis. This would provide a second very significant 
variable to use as the objective function and would 
certainly be in line with the current philosophy of "design 
Bo cost." 

The second feature would be a topsides arrangement 
routine. Naval surface combatant design has become sensi- 
tive to this area and in some cases the designs are driven 
by the requirements for deck area necessary to operate 
electronics and weapons systems. A feature that could be 
coupled with the above routine and used during the 
concept design phase would be a graphics module. There 
currently exist programs capable of generating rough hull 
forms with less information than the synthesis model is 
capable of providing. The British have developed software 
Which will interactively, with the designer, develop the 


ship's arrangements and deckplans [Ref. 24]. 
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A final possibility in the graphics area would be to 
develop a data base of three to five surface combatant hull 
forms which may be expanded, contracted, or mixed to achieve 
a desired hull form. This is similar to the approach taken 
with hydrofoil design synthesis and would provide the 
designer with additional design flexibility [Ref. 25]. 

Minor changes to the synthesis program could be made. 
Payload items used on naval combatants are continually 
changing, the items in the payload shopping list [Ref. 3], 
Should be reviewed and updated as part of a periodic program 
maintenance action. 

In conclusion, it is felt that the optimizer/synthesis 
system is a design tool which can, when properly used, 
Significantly improve the quality of feasibility studies 
without increasing the time required to accomplish the 
studies. As computer hardware and software improve, the 
design procedure outlined in this thesis will become more 
commonplace and as a result the designers will further 
benefit by designs which are developed computationally to a 
much greater degree of detail than is currently observed at 


any given phase of design. 
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APPENDIX B 


GLOBAL CATALOG 


FORTRAN DEFINITION LOCATION NAME 
DISPFL Full Load Displacement, tons 
Vsus Maximum Sustained Speed, knots 
Vend Endurance Speed, Knots 

ENDUR Endurance Range, naut. mi. 

LBP Length Btwn Perpendiculars, ft. 
ley Length to Beam Ratio 

B/H Beam to Draft Ratio 

Cp Prismatic Coefficient 

Ox Midship Section Coefficient 

G1 DISPRE/DPTRY -— 1.0 

G2 CGEEE RGR -—91.,0 

G3 2207 > Lhe — 1.0 

G4 BYR/4-0 = 1.0 

G5 -R/10.0 

G6 0243 /Cpr— 71.0 

Gi Co7O570 — 1.0 

G8 QO.5/SLRAT - 1.0 

G9 D2001,/ Cv — 150 


J0 





APPENDIX B 


GLOBAL CATALOG 


GLOBAL FORTRAN DEFINITION LOCATION NAME 
86 GEO Cv/0.006 - 1.0 

87 G11 SLRAT/2.0 - 1.0 

88 Giz -1.0 

89 G13 RSSV/DHV - 1.0 
2901 


Safer stcments 2-/6 and YO-2Zo01 are contained in the S vector 


described in Chapter II. Only those of interest to the thesis 
have been listed. 
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et 


BM 
By 
Ca 


CASDAC 


Cb 
CGFLD 
Cp 

er 

Cy 
Cwp 
Ox 
DISPFL 
DHV 
DPTRY 
FO 
BAO 
F20 
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APPENDIX C 


NOMENCLATURE 


Move Parameter 

Beam to Draft Ratio 

Beam at Midship Section, ft. 
Metacentric Radius, ft. 

Beam to Draft Ratio 

Transverse Moment of Inertia Coefficient 
Computer-aided Ship Design and Construc- 
tion 

block Coeciticrent 

VCG @ Full Load Displacement, ft. 
Prismatic Coefficient 

Residual Resistance Coefficient 
Volumetric Coefficient 

Waterplane Coefficient 

Midship Section Coefficient 

Full Load Displacement, long tons 
Deck House Volume, ft°* 
Displacement Estimate, long tons 
PEeenboardsatreouvation O, ft. 
Freeboard at Station 10, ft. 
Freeboard at Station 20, ft. 


Peco woumiace Correetion Factor, it. 
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GC ) 
GM 
GM/B 
H 

KB 

KG 
mGIRY 
LBP 
iby B 


OBJ 


RSSV 


[~” 


SC) 
SHP 
SLRAT 
VCG 


Vsus 


APPENDIX C 


NOMENCLATURE 


Constraint Values 

Metacentric Height, ft. 

Ratio of Metacentric Height to Beam 
Dit@emnoad Pratt, ft. 

Herchtwonm center Of Bouyancy, it. 
Herchte@: secentcer Of Gravity, it. 

KGa ctimacre, 1 t. 

Length Between Pervendiculars, ft. 
Length to Beam Ratio 

Objective Function 

Iteration Number 

Required Suverstructure Volume, ft° 
pearch Direction 

Meray OC OLOrIne Ship Specifications 
Shaft Horsepower 

Speed-Length Ratio 

Vertueai Center of Gravity, ft. 
Maximum Continuous Sustained Speed, kts. 


Design Variable Vector 
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